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ABSTRACT 


Pool  boiling  heat-transfer  measurements  were  made  using 
seven  Gewa-T  copper  tubes  in  R-113.   The  first  tube  (19-mm 
OD  x  0.75  fins/mm  x  0.25  mm  gap)  was  tested  under  three  condi- 
tions:  (a)  plain,  (b)  with  four  shrouds  and  (c)  with  up  to 
5  wires  (0.13  mm  diameter)  wrapped  in  the  inter-fin  trough 
area.   Verifying  prior  data,  the  shroud  with  longitudinal 
openings  of  60°  at  the  top  and  8.5°  at  the  bottom  gave  the 
best  performance  among  the  four  shrouds  tested.   This  shroud 
increased  the  boiling  heat-transfer  coefficient  by  253  per- 
cent (over  the  smooth-tube  case)  at  a  heat  flux  of  10,000 

2 
W/m  while  it  showed  only  a  18  percent  increase  at  100,000 

2 
W/m  .   When  wire  wraps  were  provided,  in  all  cases,  the 

heat- transfer  coefficient  was  improved  at  all  heat  fluxes. 
The  use  of  three  wires  gave  the  best  performance  with  341- 
percent  increase  in  heat-transfer  coefficient  at  10,000 

2  2 

W/m   and  a  130  percent  increase  at  100,000  W/m  . 

The  next  three  Gewa-T  tubes  had  dimensions  of  25  mm 

OD  x0.75    fins/mm  with    fin-tip   gaps    of    0.15,    0.25    and    0.35   mm, 

respectively.   The  last  three  Gewa-T  tubes  had  the  same 

dimensions  except  that  the  fin  density  was  1.02  fins/mm. 

For  each  of  these  fin  densities,  the  0.25  mm  fin-tip  gap 

produced  the  best  boiling  performance  at  all  heat  fluxes. 

Also,  for  a  given  fin-tip  gap,  the  boiling  performance 


increased  with  increase  in  fin  density.   For  a  fin-tip  gap 
of  0.25  mm,  the  lower  fin-density  tube  produced  an  80-percent 
increase,  while  the  higher  fin  density  tube  produced  a  103- 
percent  increase  in  the  boiling  heat-transfer  coefficient. 
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I.   INTRODUCTION 

A.   BACKGROUND 

The  energy  crisis  that  began  in  1974  has  made  a  consider- 
able impact  on  the  design  and  operation  of  energy-conversion 
devices.   Today,  designers  are  challenged  to  design  more 
and  more  efficient  heat-transfer  systems.   Among  the  numer- 
ous areas  that  designers  are  faced  with  this  challenge, 
boilers  or  evaporators  used  in  many  engineering  applications 
demand  considerable  attention.   The  use  of  advanced  boiling 
surfaces  may  lead  to  significant  improvements  in  efficiency 
and/or  significant  savings  in  the  initial  cost  and  equipment 
sizes.   Also,  due  to  the  ever-decreasing  sizes  of  electronic 
components,  the  presence  of  very  high  heat  fluxes  has  made 
the  cooling  of  such  components  a  challenging  problem.   Marto 
and  Lepere  [Ref.  1]  stated  that  the  use  of  advanced  boiling 
surfaces  may  have  potential  promise  in  the  field  of  electronic 
cooling.   However,  the  lack  of  generalized  performance  data 
and  useful  theoretical  models  for  such  advanced  surfaces 
have  still  kept  designers  away  from  the  utmost  use  of  these 
surfaces . 

Due  to  a  lack  of  complete  understanding,  the  theoretical 
treatment  of  the  nucleate  pool-boiling  process,  even  on  a 
smooth  surface  ,  is  extremely  complex.   The  existence  of  a 
large  number  of  semi-empirical  formulas  is  probably  a  very 
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good  indication  of  the  complex  nature  of  the  nucleate  boil- 
ing process.   Chongrungreon  and  Sauer  [Ref.  2]  summarized 
nine(9)  correlations  for  nucleate  pool  boiling  of  single- 
component  fluids  from  smooth  tubes.   They  stated  the  inade- 
quacy of  any  of  these  correlations  when  compared  with 
experimental  data. 

Nucleate  pool  boiling  is  strongly  dependent  on  the  sur- 
face structure;  i.e.,  the  heat  flux  increases  with  increase 
in  stable  nucleation  site  density.   The  size,  shape  and  den- 
sity of  such  nucleation  sites  are  almost  impossible  to  pre- 
dict.  However,  the  entire  nucleate  pool-boiling  process  is 
dependent  on  the  performance  at  these  nucleation  sites.   It 
is  well  known  that  the  ability  of  these  sites  to  develop  and 
trap  small  vapor  bubbles  is  the  reason  for  improved  performance 
in  nucleate  pool  boiling  in  comparison  with  natural-convective 
boiling. 

When  a  bubble  is  being  generated  at  a  surface  nucleation 
site  (or  a  pore) ,  the  wall  superheat  is  given  by  the  follow- 
ing approximate  theoretical  equation  [Ref.  3]: 


2  o  T 
T„  -  T, 


pV  hFG  R 


Thus,  for  a  given  fluid,  the  wall  superheat  can  be  de- 
creased (boiler  performance  increases  with  decrease  in  wall 
superheat)  with  increase  in  pore  radius.   Of  course,  this 
radius  must  have  a  practical  upper  limit  since  too  large  a 
pore  would  not  trap  a  vapor  bubble. 


14 


Nishikawa  and  Ito  [Ref.  4]  recently  discussed  two  methods 
to  augment  nucleate  boiling:   (1)  use  teflon-coated  pits  to 
reduce  wettability  of  boiling  liqud,  and  (2)  manufacture 
surfaces  with  numerous  re-entrant  cavities  that  have  the 
ability  to  trap  vapor  and  keep  the  nucleation  sites  active. 
Due  to  the  high  wettability  of  low  surface-tension  fluids 
(such  as  f luorcarbons) ,  the  use  of  teflon-coated  pits  has  not 
been  successful.   Thus,  the  surfaces  with  re-entrant- type 
cavities  are  the  viable  candidates  for  low  surface-tension 
liquids . 

Taking  advantage  of  the  concept  of  re-entrant  cavities, 
a  number  of  enhanced  boiling  surfaces  have  been  invented 
during  the  last  two  decades.   They  include  the  Linde  High 
Flux  surface  of  Union  Carbide,  Gewa-T  surface  of  Wieland 
and  Thermoexcel-E  surface  of  Hitachi.   All  of  these  surfaces 
are  provided  with  re-entrant  cavities  through  specialized 
manufacturing  processes.   These  re-entrant  cavities  introduce 
yet  more  parameters  to  be  desired. 

In  view  of  the  inadequate  understanding  of  the  nucleate 
pool  boiling  process  on  smooth  tubes,  a  thorough  understand- 
ing of  the  boiling  process  on  advanced  surfaces  may  take 
many  more  years  to  come.   Thus,  pbesent-day  designers  have 
very  little  information  when  considering  advanced  boiling 
surfaces.   It  is  very  important  that  researchers  gather  a 
vast  amount  of  experimental  data  covering  a  wide  range  of 
heat-transfer  fluids.   Such  extensive  data  will  enable  the 
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development  of  generalized  correlations  or  successful 
theoretical  models. 

Yilmaz,  et  al.  [Ref.  5]  compared  the  nucleate  pool-boiling 
heat- transfer  performance  of  three  copper  tubes:   Linde  High 
Flux,  Wieland  Gewa-T  and  Hitachi  Thermoexcel-E .   Their  data 
taken  with  p-xylene  showed  that  while  all  tubes  performed 
better  than  smooth  tubes,  the  best  performance  was  with  the 
Linde  High  Flux  surface.   The  Thermoexcel-E  tube  performed 
better  than  the  Gewa-T  tube  at  low  heat  flux  and  both  per- 
formed similarly  at  high  heat  flux.   These  data  were  taken 
under  the  fully-established  boiling  regime;  thus,  no  infor- 
mation was  provided  on  the  hysteresis  behavior  or  temperature 
overshoot  of  these  surfaces. 

Bergles  and  Chyu  [Ref.  6]  compared  the  nucleate  pool- 
boiling  heat-transfer  performance  of  four  different  Linde 
High  Flux  tubes  to  a  plain  copper  tube  in  distilled  water  and 
in  R-113.   They  used  three  different  methods  of  aging  the  test 
surface  prior  to  collecting  data.   They  reported  significant 
temperature  overhsoots  of  both  plain  and  High  Flux  surfaces 
in  R-113.   They  also  reported  that  these  overshoots  were 
sensitive  to  aging,  initial  subcooling  and  rate  of  power 
increase. 

Marto  and  Lepere  [Ref.  1]  obtained  results  on  the  heat- 
transfer  performance  of  three  heat-transfer  surfaces:   Linde 
High  Flux,  Thermoexcel-E,  Gewa-T,  and  a  plain  copper  sur- 
face in  R-113  and  FC-72.   They  reported  a  two  to  tenfold 
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increase  in  the  heat-transfer  coefficient  when  compared  to 
plain  surface.   The  High  Flux  surface  was  most  effective  over 
a  broad  range  of  heat  fluxes.   The  Thermoexcel-E  surface 

showed  similar  gains  in  heat-transfer  coefficient  to  that 

2 
of  the  High  Flux  surface  below  10  kW/m  .   The  Gewa-T  surface 

was  not  as  effective  as  the  other  surfaces  at  low  heat  fluxes 

2 
but  performed  comparably  at  high  heat  fluxes  (near  100  kW/m  ) 

Marto  and  Hernandez  [Ref.  7]  studied  the  boiling  per- 
formance of  a  Gewa-T  finned  surface  in  Freon-113.   They 
examined  the  influence  of  placing  aluminum  shrouds  around  the 
test  tubes,  and  found  that  the  thermal  performance  of  these 
surfaces  were  dependent  upon  the  liquid-vapor  flow  within 
the  channels  between  neighboring  T-shaped  fins.   They  also 
examined  the  influence  of  the  fins  by  progressively  machining 
away  the  fin  height  to  arrive  at  a  smooth  cylindrical  surface. 
They  reported  that  the  addition  of  the  "T-caps"  to  straight 
fins  produced  the  most  significant  improvement  in  heat- 
transfer  performance  when  compared  to  a  smooth  tube.   The 
use  of  shrouds  increased  the  heat-transfer  coefficient  at 
all  practical  heat  fluxes,  but  decreased  the  heat- transfer 
coefficient  at  very  high  heat  fluxes  when  compared  to  the 
unshrouded  surface. 

B.   THESIS  OBJECTIVES 

Based  on  the  foregoing  discussion,  the  major  motivation 
in  this  work  was  to  further  study  the  boiling  heat-transfer 
performance  of  Gewa-T  surfaces  in  R-113,  and  to  attempt 
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to  improve  this  performance  through  the  use  of  further 
enhancing  techniques. 

There  were  four  major  objectives  in  this  thesis  study: 

1.  Take  data  to  verify  results  of  Hernandez  [Ref.  8]  on 

a  standard  21.2  mm  O.D.  Gewa-T  tube  with  and  without  shrouds 

2.  Use  wire  wraps  in-  the  inter-fin  area  to  improve 
nucleate-boiling  heat-transfer  performance  from  Gewa-T 
tubes . 

3.  Take  data  using  the  combination  of  best  shroud  (found 
in  1)  and  the  best  number  of  wire  wraps  (found  in  2)  on 
Gewa-T  tubes. 

4.  Take  data  to  study  the  effect  of  fin-tip  spacing  and 
fin  density  on  the  heat-transfer  performance  of  Gewa-T 
tubes . 
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II.   EXPERIMENTAL  APPARATUS 

A.  REQUIREMENTS  GOVERNING  DESIGN 

The  purpose  of  this  research  was  the  observation  of 
nucleate  pool  boiling  from  different  Gewa-T  surfaces.   In 
order  to  get  a  more  comprehensive  knowledge  about  the  physi- 
cal mechanisms  that  govern  this  phenomenon,  it  is  necessary 
to  measure  the  parameters  involved  which  affect  the  nucleate 
boiling  process.   The  selected  parameters  to  be  obtained 
are  the  following: 

1.  Barometric  pressure 

2.  Test  section  wall  temperature 

3.  Fluid  bulk  temperature 

4 .  Vapor  temperature 

5.  Surface  heat  flux 

B.  TEST  APPARATUS 

The  experimental  apparatus  designed  by  Lepere  [Ref.  9], 
and  modified  by  Hernandez  [Ref.  8],  with  further  minor  modi- 
fications, was  used  during  this  thesis  work.   To  gather  data 
more  conveniently  at  a  faster  rate,  a  Hewlett  Packard  model 
9826  computer  and  an  HP-3054A  automatic  data-acquisition 
system  were  added  to  the  experimental  facility.   Lepere 
[Ref.  9]  and  Hernandez  [Ref.  8]  provide  detailed  descriptions 
of  the  apparatus;  therefore,  only  a  brief  description  is 
presented  in  this  thesis. 
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Figure  1  shows  a  schematic  of  the  test  apparatus.   The 
boiling  surfaces  consist  of  a  cylindrical  copper  tube  with 
a  smooth  interior  surface  and  a  commercially-available  T- 
shaped  finned  surface  manufactured  by  Wieland-Werke ,  AG.   Each 
tube  was  conditioned  for  the  insertion  of  thermocouples  by 
drilling  eight  longitudinal  holes  as  indicated  in  Figure  2 
and  introducing  the  respective  copper-constantan  thermocouples 
centered  by  means  of  capillary  copper  tubes.   In  the  center 
of  the  tube  a  cartridge  heater  was  placed  to  act  as  the  heat 
source  as  shown  in  Figure  3. 

Both  ends  of  the  test  section  were  machined  to  remove 
the  enhanced  surface  on  the  outside  and  to  increase  the 
interior  diameter.   This  reduced  the  cross  sectional  area 
and  thereby  minimized  the  heat  transfer  from  the  ends  of  the 
test  section.   A  Teflon  plug  was  inserted  into  the  free  end 
for  sealing  and  heat  transfer  reduction  purposes.   The  other 
end  of  the  test  section  was  connected  to  a  long  copper  ex- 
tension tube  through  a  90-degree  elbow,  which  allowed  the  test 
section  to  be  placed  in  a  horizontal  position  within  a  thick- 
walled  Pyrex  glass  vessel.   Since  one  of  the  requirements 
of  the  apparatus  was  to  maintain  the  liquid  at  saturation 
temperature  and  atmospheric  pressure,  a  plate  heater  and  a 
vent  tube  were  provided  to  satisfy  these  conditions.   In 
addition  avoiding  excessive  loss  of  vapor,  a  water-cooled 
copper  coil  was  placed  in  the  vapor  space  of  the  vessel,  to 
act  as  a  primary  condenser.   The  vessel  was  covered  with  a 
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Figure    1.       Schematic    of   Test   Apparatus 
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Plexiglas  lid,  which  had  been  fitted  with  a  rubber  O-ring. 
Several  holes  were  drilled  in  the  lid  to  fix  in  place  the 
test  section  extension,  vent  tube,  connections  for  primary 
and  secondary  condensers,  and  two  more  thermocouples  to  sense 
both  the  bulk  liquid  and  vapor  temperatures.   The  Pyrex 
glass  vessel  with  the  attached  fittings  was  placed  on  the 
plate  heater,  which  served  to  maintain  the  liquid  at  satura- 
tion temperature.   A  secondary  condenser  was  mounted  external 
to  the  vessel. 

C.  INSTRUMENTATION 

The  power  supply  was  connected  to  the  cartridge  heater 
through  a  voltage  regulator  and  then  to  a  Variac  scaled  from 
0  to  120  VAC.   This  arrangement  allowed  fine  adjustment  for 
a  stable  feed  voltage  during  the  runs.   A  precision  resistor 
of  2.031  ohms  was  connected  in  series  with  the  cartridge 
heater  to  provide  a  measure  of  the  current  flow  by  reading 
the  voltage  drop  across  this  element.   A  digital  voltmeter 
was  used  to  measure  both  the  voltage  across  the  resistor  and 
the  voltage  across  the  cartridge  heater.   All  the  thermo- 
couples were  connected  to  a  computer  data  acquisition  system, 
and  ambient  conditions  and  voltage  readings  were  supplied 
manually.   The  microcomputer  processed  all  this  information 
through  a  program  included  in  Appendix  A. 

D.  CALIBRATION  OF  THERMOCOUPLES 

A  total  of  fifteen  wall  thermocouples  in  tubes  1  and  7 
(one  thermocouple  in  tube  1  was  defective,  hence  it  was 
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excluded) ,  and  the  two  thermocouples  for  vapor  and  liquid 
temperature  measurements  were  calibrated  in  an  insulated, 
stainless-steel  bath.   The  tubes  were  completely  immersed 
into  the  bath,  which  was  provided  with  a  motor-driven  mixer. 
The  bath  temperature  was  measured  using  a  platinum-resistance 
thermometer  accurate  to  ±0.01  K. 

Since  the  temperature  measurements  during  this  investiga- 
tion ranged  from  45°C  to  60 °C,  thermocouple  calibration  was 
performed  in  the  range  of  45°C  to  70°C.   To  begin  the  cali- 
bration process,  hot  water  (at  about  75°C)  was  added  to  the 
insulated  container.   Following  proper  mixing,  all  thermocouple 
readings  were  automatically  recorded  through  the  data- 
acquisition  system.   The  bath  temperature  was  measured  using 
the  platinum- resistance  thermometer.   To  obtain  other  lower 
temperatures,  small  quantities  of  water  were  added  to  the 
bath,  allowing  5  minutes  for  temperature  stabilization  of  the 
test  section. 

The  discrepancy  (i.e.,  the  platinum- resistance  thermometer 
reading  minus  the  thermocouple  reading)  was  plotted  against 
the  thermocouple  reading  as  shown  in  Figure  4 .   The  second 
order  polynomial  curve  shown  in  this  figure  was  generated 
using  the  data  from  all  seventeen  thermocouples.   However, 
for  clarity,  only  three  thermocouples  were  selected  for  data 
points  from  each  tube.   It  is  evident  that  the  thermocouples 
in  both  tubes  show  similar  discrepancies,  the  thermocouple 
readings  have  only  about  ±0.1  K  discrepancy  around  the 
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calibration  line.   Based  on  this  fact,  no  calibration  was 
performed  for  the  other  five  tubes.   Instead,  the  calibration 
curve  generated  for  tubes  1  and  7  was  assumed  for  all  tubes. 
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Ill .   EXPERIMENTAL  PROCEDURES 

A.   NORMAL  OPERATION 

The  Procedure  B  employed  by  Marto  and  Hernandez  [Ref.  7] 
has  been  demonstrated  to  be  the  most  adequate  to  carry  out 
the  boiling  runs.   Because  of  the  large  influence  of  surface 
past  history  on  boiling  incipience,  this  procedure  was  estab- 
lished to  allow  repeatability  and  permit  comparison  of  the 
obtained  data  for  different  enhancements. 

The  procedure  is  as  follows : 

1.  Once  the  test  section  has  been  immersed  in  the  pool, 

it  was  subjected  to  one  hour  of  pre-boiling  by  setting  the 

2 
cartridge  heater  to  give  a  heat  flux  of  about  30  kW/m  .   The 

plate  heater  was  adjusted  to  the  maximum  temperature  for  this 

one-hour  period.   This  initial  vigorous  boiling  served  to 

degas  the  liquid  and  to  force  the  noncondensable  gases  out  to 

the  atmosphere. 

2.  After  the  initial  aging  process,  both  test-section  heater 
and  plate  heater  were  secured,  allowing  cooling  of  the  liquid 
and  the  test  section  for  about  30  minutes. 

3.  The  plate-heater  voltage  was  then  adjusted  in  order  to 
maintain  the  liquid  at  saturated  conditions  and  the  cartridge 
heater  was  turned  on  to  the  first  setting  of  12  volts. 

4.  For  all  the  consecutive  settings  during  both  increasing 
and  decreasing  heat  fluxes,  the  boiling  was  allowed  to 
stabilize  for  five  minutes  at  each  power  setting. 
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5.   At  a  given  power  setting,  the  following  data  were 
recorded:   heater  voltage,  precision-resistor  voltage,  vapor 
temperature,  liquid  bulk  temperature  and  wall  temperatures 
of  the  test  section. 

B.   GEWA-T  SURFACES 

Three  groups  of  surfaces  were  tested.   The  first  surface 
was  the  same  as  used  by  Hernandez  [Ref.  8]  with  a  fin  density 
of  0.75  fins/mm,  an  outside  diameter  of  21.2  mm,  and  an 
inter-fin  gap  of  0.25  mm.   The  second  group  of  three  tubes 
had  0.75  fins/mm,  an  outside  diameter  of  23.1  mm  and  inter- 
fin  gaps  of  0.15,  0.25  and  0.35,  and  the  third  group  of  three 
tubes  had  1.02  fins/mm  and  all  other  specifications  were  the 
same  as  the  second  group. 

The  geometry  of  the  enhancement  of  these  seven  Gewa-T 
surfaces  was  observed  by  making  a  longitudinal  cut  of  the 
tubes  and  photographing  them  with  50X  magnification  (Figure  5] 
As  can  be  seen,  the  tubes  with  0.75  fins/mm  had  an  approxi- 
mately circular  cavity  shape,  while  the  1.02  fins/mm  tubes 
had  a  somewhat  elongated  cavity  shape.   This  difference  in 
shape  was  probably  due  to  the  technique  used  in  construction 
to  obtain  the  various  specifications.   All  the  T-shaped  fins 
had  a  fin  height  of  approximately  0.9  5  mm. 

The  exact  determination  of  diameters  at  the  base  of  the 
fins  was  carried  out  from  the  magnified  photographs  of  the 
cavities.   This  diameter  was  determined  by  subtracting  twice 
the  fin  height  from  the  measured  outside  diameter: 
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(a)   Gewa-T  Tube  Number  2 


(b)   Gewa-T  Tube  Number  5 


Figure  5.   Optical  Micrograph  of  Gewa-T  Fin  Profile,  50X 
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D2   =   O.D.  -  2HF  , 


resulting  in  the  dimensions  presented  in  Table  1 

TABLE  1 
Dimensions  of  Gewa-T  Finned  Tubes 


UBE 
MO. 

FINS/ 
mm 

FIN  TIP 
GAP  (mm) 

FIN  HEIGHT 
(mm) 

O.D. 
(mm) 

D2 
(mm) 

1 

0.75 

0.25 

0.865 

21.2 

19.47 

2 

0.75 

0.15 

0.980 

24.99 

23.05 

3 

0.75 

0.25 

0.960 

24.99 

23.07 

4 

0.75 

0.35 

0.970 

24.99 

23.05 

5 

1.02 

0.15 

1.050 

25.30 

23.20 

6 

1.02 

0.25 

1.050 

25.30 

23.20 

7 

1.02 

0.35 

1.125 

25.30 

23.05 

For  simplicity,  the  diameter  at  the  base  (D~)  was  assumed 
to  be  19.5  mm  for  Tube  1  and  23.1  mm  for  all  other  tubes. 

C.   HEAT-FLUX  CALCULATION 

According  to  earlier  description  about  test-section  con- 
struction, the  cartridge  heater  acts  as  the  heat  source, 
inserted  into  the  tube,  and  variations  in  heat-transfer  rate 
are  made  through  different  power  settings  with  a  variac. 


QH  =  vi 
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1) 


where : 


Q    =   heat-transfer  rate  from  the  cartridge 
heater  (W) 

V   =   voltage  across  the  cartridge  heater 
element  (volts) 

I    =   current  through  the  heater  element  (amps 


This  circuit  is  connected  in  series  with  a  precision 
resistor,  which  provides  the  means  to  compute  the  electrical 
current  by  the  measurement  of  the  voltage  drop  across  the 
resistor : 


1   =   VRES/R  (2) 


where : 

vt3tto   =   voltage  drop  across  the  precision 
resistor  (volts) 

R     =   resistance  of  the  precision  resistor  (ohms) 

=   2.0  31  ohms 

Both  ends  of  the  Gewa-T  tube  were  machined  to  provide  a 
smooth  surface  and  a  reduced  wall  thickness  in  order  to  mini- 
mize longitudinal  heat  conduction  and  therefore  the  heat 
losses  from  the  unenhanced  ends.   The  geometry  of  these  ends 
is  shown  in  Figure  6,  where: 

b    =   base  of  the  T-shaped  fins 

T,   =   temperature  at  the  base  of  T-shaped  fins 
13      (°C) 
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Figure  6.   Geometry  of  Unenhanced  Ends 


For  this  cylindrical  geometry: 


Ac  -  \<-^ 


(3) 


Ac   =   Px  L 
S         c 


(4; 


where 


=   cross-sectional  area  of  tube  wall  (m  ) 


D    =   tube  outside  diameter  (m) 


=   tube  inside  diameter  (m) 


A„   =   tube  outside  surface  area  (m' 


P    =   tube  outside  wall  perimeter  (ml 


It  was  assumed  that  the  temperature  at  the  base  of  the 
thin-walled  fin  was  equal  to  the  temperature  at  the  base  of 
the  T-shaped  fins.   Therefore: 
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(x)   =   T(x)  -  T 


SAT 


5) 


In  general , 


d  9     2 

— j    -   m 

dx^ 


=   0 


where : 


m 


h  P 
kA^ 


7) 


Therefore 


QF   =   QB 


i  *   dT 
C  dx 


=   -kA 


x=0 


de 

C  dx 


x=0 


Now,  assuming  negligible  convective  heat  loss  from  the 
thin  tube  tip: 


de_ 

dx 


=      0 


x= 


^ 


(9 


and 


0  Cosh    m(L-x) 


B 


Cosh   inL 


(10 


For    this    case 


-      fi 


PkcAc     9B    Tanh    ml^ 


11) 
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'B   =   S<°»   ■   TB  -  TSAT  <12 


where 


Q     =   heat-transfer  rate  through  unenhanced 
end  (W) 

TqAT   =   fluid  saturation  temperature  (°C) 

k      =   thermal  conductivity  of  the  copper  (W/m.K) 

_  2 

h      =   average  heat-transfer  coefficient  (W/m  .K) 


Assuming  the  wall  temperature  is  the  same  as  the  average 
of  the  eight  calibrated  wall  temperature  thermocouple 
measurements,  less  temperature  drop  due  to  radial  conduction 
from  thermocouple  position  to  the  base  of  the  T-shaped  fin, 


ln(D2/D1) 


TB   =   TAVG  "  QH   2,Lk  <13> 


T      =   (  I      T  )/8.0  (14) 

AVG       n41   n 


where 


D,   =   diameter  at  the  position  of  thermocouples 
1      (m) 

T  vr   =   average  wall  temperature  (°C) 
For  the  geometry  shown  in  Figure  6 : 


Lc  .=   Lu  +  |  (15) 
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t   =   D   -  D.  (16) 

o     1 


where 


L    =   corrected  length  (m) 

L    =   length  of  unenhanced  ends  (m) 

t    =   tube  thickness  (m) 

The  average  difference  between  the  wall  temperature  and 
the  saturation  temperature  may  now  be  determined  from  the 
following  equation: 


..  C  Cosh  m(Lr-x) 

=   T="   /     ~n u T^ dx  (17) 

Lc  ^     Cosh  mLc 


Tanh  mLp  (18) 


mL„         C 


"   TW  -  TSAT  (19 


where 


(T  -T    )   =   average  difference  between  surface 
wall  temperature  and  fluid 
saturation  temperature  (K) 

Since  the  unenhanced  ends  are  very  long  and  the  thickness 
is  very  small,  only  very  little  heat  transfer  would  take 
place  at  these  ends.   Thus,  these  ends  undergo  free  convection, 
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for  which  the  Churchill  and  Chu  correlation  as  stated  in 
[Ref.  10]  for  the  average  Nusselt  number  was  assumed: 


_  1/6 

j  0.3  87  RaD 

nu   =  ,0.60+ °  ; 

Do  (  [1  +  (0.559/Pr)9/16]8/27 


-5         —  12 

10    D    <    RaD      <    10 

o 


The  average  Nusselt  number  is : 


h  D 
NuD    -   -^  (2i; 

o 


where : 

k   =   thermal  conductivity  of  Freon-113  (W/m.K) 
Solving  for  h: 


0.387  Ra1/6 


I 


h    =    A    0.60+ °  ; 

Do  I       [1  +  (0.559/Pr)9/16]8/27 


Pr   =   ^  (23) 


2 
q  S  (T   -  T    )  D 
d.  W    SAT;  o  ,_. 

Ra     =   24 

D  v  a 


1  Ap 
p  AT 


25 


37 


where : 


Pr   =   Prandtl  number 

2 

v    =   kinematic  viscosity  (m  /s) 

2 
a    =   thermal  diffusivity  (m  /s) 

Ra     =   average  Rayleigh  number 


g 

b 


gravitational  acceleration  (m/s  ) 

volumetric  thermal  expansion  coefficient 
(1/K) 


p    =   density  of  freon-113  (kg/m  ) 

Now,  substitution  of  Equations  7,  18  and  24  into  Equation  22 
results  in: 


0.387 


hP  ^2   u  1/6 


gBDoeBTanh{(j^-)    Lc 


vaL 


hP  » 1/2 


(26) 


C  vkA, 


h   = 


D   \ 
o 


0.60  + 


[1  +  (0.559/Pr)9/16]8/27 


Equation  26  is  solved  for  h  by  iterative  technique 
within  a  range  of  precision  of  0.001.   Knowing  the  value  of 
natural-convective  heat-transfer  coefficient  at  unenhanced 
ends,  the  total  heat-loss  rate  is  calculated  from  Equation 
11,  and  the  heat-transfer  rate  through  the  enhanced  Gewa-T 
surface  is  obtained  by  subtracting  the  total  heat-loss  rate 
from  the  total  heat-transfer  rate: 


QLOSS   =   2QF 


(27) 
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Q     =     QE~   2QF 


28 


Finally,  the  heat  flux  from  the  Gewa-T  surface  at  the 
base  of  the  T-shaped  fins  is: 


Q/^ 


(29) 


Ag   =   tt  D„  L 


(30) 
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IV.   RESULTS  AND  DISCUSSION 

A.   DESIGNATION  OF  DATA  FILES 

A  total  of  seven  Gewa-T  tubes  were  tested  during  this 
investigation.   As  can  be  seen  from  Table  1,  the  primary 
geometric  variables  were:   fin-tip  gap,  fin  density  and 
outside  diameter  of  tubes  including  fins. 

The  general  designation  for  data  files  is: 

GTNPnn 

where : 

GT   indicates  Gewa-T  surface 

N    specifies  the  number  of  the  tube  tested 
(Table  1) 

P    indicates  the  position  of  the  test  section 

(i.e.,  H  for  horizontal  arid  V  for  vertical), 
and 

nn   indicates  the  corresponding  cumulative  run 
number . 

For  example,  GT5H28  specifies  tube  number  five  tested  in  the 
horizontal  position  during  run  number  28.   From  Table  1, 
the  corresponding  specifications  for  this  test  section  are: 
25.3  mm  outside  diameter,  23.1  mm  diameter  at  the  base  of  the 
fins,  1.02  fins/mm  and  inter-fin  gap  of  0,15  mm.   While  the 
file  name  just  presented  contain  data  in  the  very  raw  form, 
the  processed  data  (i.e.,  q  versus  AT)  were  stored  in  a 
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separate  file.   This  plot  data  file  could  easily  be  recognized 
by  the  "P"  added  at  the  beginning  of  the  raw  data  file  name. 
In  this  manner,  the  plot  data  file  corresponding  to  the  above 
raw  data  file  has  the  name  PGT5H28. 

A  total  of  thirty  different  data  runs  were  carried  out 
during  this  investigation.   These  runs  break  down  in  the 
following  manner: 

Run  Number        Tube  Number        Test  Condition 

Plain  -  Vertical  position 
Plain  -  Horizontal 
With  shrouds  and/or  Wires 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 
Plain  -  Horizontal 

B.   CIRCUMFERENTIAL  UNIFORMITY  OF  HEAT  FLUX  AROUND  HEATER 
ELEMENT 

During  the  early  stages  of  this  investigation,  it  was 

thought  necessary  to  test  if  the  heater  element  would  provide 

a  non-uniform  heat  flux  along  the  circumferential  direction. 

Such  a  non-uniform  heat-flux  variation  was  assumed  possible 

based  on  the  heater-element  construction.   Therefore,  tube  1 

was  immersed  vertically  in  the  R-113  liquid  pool  and  three 
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runs  were  made.   All  of  these  runs  showed  only  a  negligible 
variation  of  the  wall  thermocouple  temperatures;  maximum 
variation  was  ±0.25  K.   These  runs  clearly  showed  that  the 
heater  element  indeed  provided  a  uniform  heat  flux  along  the 
circumferential  direction.   No  detailed  data  are  presented 
in  this  thesis  for  the  tube  in  the  vertical  position.   All 
other  results  presented  herein  are  for  the  tubes  in  the 
horizontal  position. 

C.  LONGITUDINAL  UNIFORMITY  OF  HEAT  FLUX  ALONG  THE  TEST 
SECTION 

During  most  of  the  data  runs  with  the  tubes  in  the  hori- 
zontal position,  somewhat  unexplainable  phenomena  were 
observed.   As  the  heat  flux  was  gradually  increased,  bubble 
formation  began  simultaneously  at  the  two  ends,  while  the 
center  portion  of  the  tube  underwent  only  natural  convection. 
Approximately  40  percent  of  maximum  heat  flux  was  necessary 
to  obtain  nucleate  boiling  from  the  entire  tube.   Since  all 
the  wall  thermocouples  are  located  at  midway  along  the  tube, 
the  average  wall  temperature  values  provided  in  all  data 
runs  may  be  in  slight  error. 

D.  PLOT  ANALYSIS  FOR  PLAIN  SURFACE 

Figure  7  represents  a  typical  curve  for  nucleate  pool- 
boiling  process  of  Gewa-T  surface  in  Freon-113 .   The  observed 
behavior  of  this  process  is  analyzed  in  this  section  by 
studying  the  different  sections  of  the  boiling  curve. 
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Gewa-T  Surface 
in  Freon-113 


B  +■  Incipient  Boiling  Point 


AT 


A-B 
B-C 

C-D 
D-E 

E-F 


Natural  Convection 

Mixed  Boiling  (Activation  of 

Nucleation  Sites) 

Active  Boiling 

High  Flux  Nucleate  Boiling 

Stable  Nucleate  Boiling 


Figure  7.   Typical  Nucleate  Pool-Boiling  Curve 
for  a  Gewa-T  Surface 
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From  point  A  to  point  B,  a  continuous  increase  in  heat 

flux,  accompanied  by  an  increase  in  wall  superheat  (TTT  -  T_._) 

W     SAT 

is  observed  when  heat  flux  is  increased.   This  region  corres- 
ponds to  natural-convection  process,  as  also  verified  during 
experimental  runs  with  the  absence  of  bubbles  along  both 
enhanced  surface  and  unenhanced  ends . 

From  point  B  to  point  C,  a  drastic  reduction  in  wall 
superheat  was  observed  while  the  heat  flux  continuously 
increased.   This  region  was  somewhat  unstable  as  more  and 
more  nucleation  sites  became  active.   This  region  is  known 
as  the  mixed  boiling  region,  where  transition  from  natural- 
convective  boiling  to  nucleate  pool  boiling  takes  place.   At 
these  heat  fluxes,  visual  observations  revealed  that  the 
bubble  formation  was  quite  random.   While  some  portions  of 
the  tube  showed  bubbles,  other  portions  showed  no  bubble 
formation.   The  wall  superheat  continued  to  decrease  until 
all  nucleation  sites  became  active,  as  indicated  by  point  D. 

Once  all  available  nucleation  sites  have  been  activated, 
the  wall  superheat  starts  to  increase  with  increasing  heat 
flux  as  shown  by  region  D  to  E.   When  boiling  in  this  region, 
the  entire  tube  showed  very  high  density  of  bubble  formation. 
At  this  point,  it  is  important  to  point  out  that  during  all 
experimental  runs,  for  the  entire  range  of  power  settings  (7.5  to 
485  W) ,   little  or   no  nucleate  pool  boiling  on  the  unenhanced 
surface  of  the  ends  was  observed  (see  Figures  8  and  9)  for 
the  different  test  sections  tested.   This  observation 
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confirmed  the  validity  of  the  mathematical  model  assumed  for 
data  reduction.   However,  a  small  rate  of  bubble  formation 
was  evident  at  the  tip  of  the  unenhanced  fin  end;  these 
bubbles  were  in  fact  generated  at  the  base  of  the  fin  (i.e., 
between  the  Teflon  plug  and  the  inner  surface  of  the  fin) . 

The  vapor  bubbles  generated  in  the  inter-fin  cavities 
are  carried  up  into  the  channels  around  the  tube  due  to 
buoyancy  forces.   As  these  bubbles  travel  up,  they  activate 
more  and  more  nucleation  sites  as  well  as  they  help  for  the 
removal  of  larger  bubbles  from  the  channels.   This  phenomenon 
is  generally  referred  to  as  the  "Chimney  Effect." 

If  the  heat  flux  is  too  large  (see  Figure  7) ,  the  rate  of 
vapor  bubble  generation  may  be  too  large  that  the  interfin 
channels  could  be  full  of  a  vapor  blanket  that  may  disrupt 
the  flow  of  liquid  to  the  hot  metal  surface  where  active 
boiling  takes  places.   At  this  point,  further  increase  in 
heat  flux  may  result  in  excessive  wall  temperatures.   Under 
these  conditions,  vapor  bubbles  were  seen  ejecting  radially 
outward  from  the  test  surface.   Line  D-E  represents  the  boil- 
ing regime  just  described.   It  is  clear  that  the  slope  (i.e., 
heat- transfer  coefficient)  at  E  is  considerably  smaller  than 
that  at  point  D.   This  decrease  in  heat-transfer  coefficient 
is  assumed  to  be  due  to  the  existence  of  an  insulating  vapor 
blanket. 

When  heat  flux  is  gradually  decreased,  the  curve  follows 
a  different  path  (D-F)  as  shown  in  Figure  7.   This  is  due 
to  the  existence  of  stable  nucleation  sites  that  will  remain 
active  for  a  wide  range  of  heat  fluxes. 
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E.  REPEATABILITY  WITH  PREVIOUS  DATA 

Following  modifications  made,  especially  to  the  data 
gathering  means  (i.e.,  automatic  data-acquisition  system  and 
newly-written  computer  programs) ,  it  was  necessary  to  take 
data  for  tube  1  (this  is  the  same  tube  tested  by  Hernandez) 
to  check  the  present  data  with  Hernandez . data  under  similar 
conditions.   Even  though  a  total  of  5  different  runs  were 
made  during  this  stage  through  debugging  of  programs  and 
system  familiarization,  only  run  number  5  will  be  presented 
in  this  thesis  with  a  comparison  to  Hernandez  data. 

Figure  10  shows  the  comparison  between  the  present  data 
and  the  data  of  Hernandez.   A  very  good  agreement  is  evident 
as  the  two  sets  disagree  by  less  than  2  percent  for  the 
nucleate  boiling  regime.   However,  considerable  disagreement 
exists  in  the  region  of  increasing  heat  flux  before  total 
nucleate  boiling.   This  is  believed  to  be  due  to  the  rather 
random  behavior  of  the  boiling  process  during  this  regime 
as  it  will  be  presented  in  the  next  section,  even  two  consecu- 
tive runs  for  a  given  tube  would  not  give  the  same  results 
during  the  unstable  boiling  regime. 

F.  REPEATABILITY  STUDY  OF  PRESENT  WORK 

Figure  11  shows  data  from  two  different  runs  completed 
on  two  different  days.   As  discussed  earlier,  a  considerable 
disagreement  can  be  seen  before  the  tube  undergoes  complete 
nucleate-boiling  process.   However,  for  the  complete 
nucleate-boiling  process,  the  data  from  different  data  runs 
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agree  to  within  0.2  percent  at  a  heat  flux  of  30,000  to 

2 
100,000  W/m  .   These  data  runs  show  considerably  good  repeata- 
bility of  this  experiment.   (Similar  agreements  were  also 
found  for  all  other  tubes.) 

Also,  these  data  show  a  considerable  enhancement  of 
Gewa-T  tube  performance  over  the  "Smooth"  tube  data  of 
Hernandez.   For  example,  at  AT  =  3.1  K,  the  Gewa-T  tube 
gives  a  4.65  times  larger  heat  flux  than  the  smooth  tube. 

As  also  discussed  by  Lepere  [Ref.  9]  and  Hernandez  [Ref. 
8]  this  enhancement  in  boiling  performance  is  mainly  attri- 
buted to  the  existence  of  re-entrant  cavities  in  the  Gewa-T 
tube. 

G.   BOILING  PERFORMANCE  OF  SMOOTH  TUBE 

To  compare  the  performance  of  Gewa-T  tubes  to  smooth-tube 
performance,  data  of  Hernandez  were  chosen.   It  is  important 
to  point  out  that  his  smooth  tube  did  not  consist  of  a  commer- 
cially-available surface.   Instead,  Hernandez  obtained  his 
"Smooth"  surface  by  machining  fins  off  of  a  Gewa-T  tube. 
He  did  not  remove  fins  completely;  the  remaining  portions  on 
the  tube  measured  about  0.1  mm. 

Figure  12  shows  Hernandez 's  data  (heat  flux  as  a  function 
of  wall  superheat)  for  only  decreasing  heat  fluxes.   A 
least-squares- fit  curve  generated  for  these  data  points  will 
be  used  in  future  figures  representing  data  on  tube  1.   Since 
tubes  2  through  7  have  a  larger  diameter,  it  was  necessary 
to  revise  Hernandez's  data  to  include  the  difference  in 
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diameter.   For  this  purpose,  a  semi-empirical  correlation 
developed  by  Cornwell  et  al.  [Ref.  11]  is  used: 


NUTB   =   CTB  ReTB3 


hTB  D   =         qD    2/3 
kF         TB  hFG  V 


PFDVA      mD      GD 


G   =   q/h. 


^F  A      PFA     yF  '  *'  FG 


CTB  kFfhTBATD//3 
hTB  D  -(hpGyF 


1/3      _TB_^F,AT Dx  ' 

hTB  D   (hFG  pp> 


r   k   7         ? 
,         ,    TB  KF,    ,  ATD  , 

hTB   =   (— ^-}    (hF7U;) 


hTB   "   (CTB  V    (h~~  yj     D 

r  (a   r 


Using  this  relationship,  Hernandez's  data  were  revised 
and  a  least-squares-fit  curve  was  generated.   This  curve  is 
also  plotted  in  Figure  12. 

H.   BOILING  PERFORMANCE  OF  GEWA-T  TUBE  WITH  SHROUDS 

The  same  four  shrouds  tested  by  Hernandez  were  used 
during  this  investigation.   These  shrouds  had  the  following 
geometry : 
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Shroud  Ni 

imber 

Upp< 

sr  Opening 

Lowe 

r  Opening 

1 

60° 

60° 

2 

60° 

30° 

3 

30° 

30° 

4 

60° 

8.5° 

All  of  these  shrouds  had  been  designed  for  tube  1;  thus, 
no  shrouds  were  used  for  tubes  2  through  7 . 

Figures  13,  14,  15,  and  16  show  data  on  tube  1  with  shroud 
numbers  1,  2,  3  and  4,  respectively.   As  also  reported  by 
Hernandez  [Ref.  8],  all  the  shrouds  improved  boiling  per- 
formance at  low  heat  fluxes   compared  to  plain  Gewa-T  tube. 
Table  2  lists  performance  calculations  for  tube  1  with  various 
shrouds.   It  can  be  seen  that  the  60°  *8.5°  shroud  (number  4) 

produced  the  maximum  improvement  of  253  percent  in  the  boiling 

2 
coefficient  over  smooth  tube  at  a  heat  flux  of  10,000  W/m  . 

At  this  heat  flux,  the  30°  x30°  shroud  resulted  in  a  165 

percent  improvement  in  the  boiling  heat-transfer  coefficient. 

2 

At  a  heat  flux  of  50,000  W/m  ,  all  shrouds  resulted  in  about 

190  percent  improvement,  except  the  30°  x30°  shroud  which 

resulted  in  a  180  percent  improvement. 

2 

At  very  high  heat  fluxes,  for  example  at  q  =  100,000  W/m  , 

the  60°  x60°  shroud  showed  a  maximum  performance  increase  of 
106  percent;  the  60°  x30°  shroud  resulted  in  a  49  percent 
improvement;  the  60°  *8.5°  shroud  resulted  in  a  18  percent 
increase;  and  the  30°  x30°  shroud  resulted  in  an  13  percent 
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improvement.   It  can  be  seen  that  the  values  just  presented 
are  somewhat  proportional  to  the  summation  of  the  two  shroud 
angles.   This  behavior  indicates  that  at  very  high  heat 
fluxes,  too  small  a  shroud  opening  can  obstruct  the  outflow 
of  vapor  and/or  the  inflow  of  liquid. 

However,  it  must  be  noted  that  at  all  practical  heat 
fluxes  (10,000  to  80,000  W/m  )  ,  all  the  shrouds  resulted  in 
a  considerable  enhancement  in  the  boiling  performance.   The 
60°  *8.5°  shroud  appears  to  be  the  best  selection  for  all 
practical  heat  fluxes. 

I.   PERFORMANCE  OF  GEWA-T  TUBE  WITH  WIRE  WRAPS  IN  INTER-FIN  GAP 

A  series  of  runs  were  made  with  two  to  five  0.1  mm-diameter, 
copper  wire  wraps  in  the  inter-fin  gap  of  tube  1.   Figures 
17,  18,  19  and  20  show  data  with  wire  wraps  of  2 ,  3,  4  and 
5,  respectively.   It  can  be  seen  that  all  of  these  curves 
show  increases  in  boiling  performance  at  all  wall  superheats 
for  decreasing  heat  flux  values.   This  enhancement  is  most 
probably  due  to  the  increase  in  the  nucleation  sites  within 
the  inter-fin  cavities.   However,  there  is  a  small  reduction 
in  performance  in  the  natural  convection  process  probably 
owing  to  that  wires  could  act  as  an  insulating  material  due 
to  the  thermal  contact  resistance  between  the  wires  and  the 
inter-fin  area. 

Table  2  also  shows  that  the  use  of  three  wires  resulted 
in  the  best  performance  increase.  This  combination  resulted 
in  a  341  percent  improvement  in  the  boiling  heat-transfer 
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coefficient  over  the  smooth  tube  value  at  a  wall  superheat 
of  1.4  K. 

J.   BOILING  PERFORMANCE  OF  TUBE  1  WITH  60°  x8.5°  SHROUD  AND 
THREE  WIRES 

Data  were  also  taken  to  combine  the  best  shroud  and  the 

best  number  of  wire  wraps  in  the  inter-fin  channel.   For 

this  purpose,  three  wire  wraps  were  added  to  the  inter-fin 

channel  and  data  were  taken  with  the  60°  x  8 .  5°  shroud.   As 

Figure  21  and  Table  2  show,  this  combination  resulted  in  a 

performance  yet  superior  to  the  use  of  three  wires  or  the  use 

of  60°  x8.5°  shroud  separately.   Table  2  shows  that  this 

arrangement  resulted  in  481  percent  increase  in  the  boiling 

heat-transfer  coefficient  over  the  smooth-tube  value  at  a 

heat  flux  of  10,000  W/m2 . 

K.   EFFECTS  OF  FIN-TIP  GAP  AND  FIN  DENSITY  OF  GEWA-T 
TUBES  ON  BOILING  PERFORMANCE 

As  described  earlier,  six  Gewa-T  surfaces  were  tested 
during  this  stage  of  the  investigation;  tubes  2,  3,  and  4 
had  a  fin  density  of  0.75  fins/mm,  and  tubes  5,  6  and  7 
with  a  fin  density  of  1.02  fins/mm,  and  each  set  of  tubes 
with  fin-tip  gaps  of  0.15,  0.25  and  0.33  (see  Table  1). 

Micrographs  were  obtained  for  all  enhanced  surfaces  on 
a  Scanning  Electron  Microscope  with  a  magnification  of  50X. 
Figure  5  shows  micrographs  representing  tubes  2  and  5. 
Computed  surface  areas  and  area  ratios  (with  respect  to 
smooth  tube)  are  listed  in  Table  3. 
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TABLE  3 
Computed  Surface  Areas  and  Area  Ratios 


Tube  Cavity  and  Fin2     Total  Surface  Area    AF/A< 
Surface  Area  (m  ) 

1  240.47  x 10"6 

2  288.01  x 10"6 

3  305.34  x 10"6 

4  285.07  x 10~6 

5  256.04  x 10~6 

6  248.61  x 10"6 

7  267.91  x 10"6 


(mz 

) 

8.818 

X 

10-3 



10.944 

X 

ID"3 

2.74 

11.603 

X 

ID"3 

2.91 

10.833 

X 

10-3 

2.72 

13.314 

X 

ID"3 

3.60 

12.928 

X 

10-3 

3.50 

13.931 

X 

ID"3 

3.77 

It  can  be  seen  that  tubes  with  a  higher  fin  density 
(1.02  fins/mm)  result  in  larger  area  ratios  than  the  tubes 
with  lower  fin  density  (0.75  fins/mm) .   However,  for  a  given 
fin  density,  the  variation  of  the  area  ratio  with  fin-tip 
gap  is  not  conclusive. 

Figures  22  through  27  show  data  taken  on  tubes  2  through 
7,  respectively.   These  results  are  also  listed  in  Table  2. 
This  table  shows  that  tube  6  (with  a  fin  density  of  1.02 
fins/mm  and  a  fin-tip  gap  of  0.25  mm)  resulted  in  the  highest 
heat  flux  at  all  values  of  wall  superheat. 

To  investigate  the  effects  of  fin  density  and  fin-tip 
gap  on  boiling  heat-transfer  performance,  the  results  pre- 
sented in  Figures  22  through  27  (which  showed  variations  of 
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heat  flux  as  a  function  of  wall  superheat)  were  re-plotted 
on  a  different  basis  as  shown  in  Figure  28.   This  figure  shows 
the  variations  of  the  enhancement  of  heat-transfer  coefficient, 
over  the  smooth-tube  value,  with  wall  superheat  for  all  six 
tubes.   It  is  evident  that  tube  6,  which  has  a  fin  density 
of  1.02  fins/mm  and  a  fin-tip  gap  of  0.25  mm  results  in  the 
maximum  enhancement  at  all  wall  superheat  values. 

Tube  5  shows  poorer  performance  than  tube  6.   This  poorer 
performance  is  believed  to  be  due  to  the  fact  that  the  smaller 
fin-tip  gap  disrupts  the  flow  of  liquid  into  the  channels 
and/or  the  vapor  flow  out  of  the  channels.   Tube  7  (which  has 
the  largest  fin-tip  gap  of  0.35  mm)  shows  the  poorest  per- 
formance among  the  set  of  these  tubes  with  a  fin  density  of 
1.02  fins/mm.   As  visually  observed  during  experimental 
runs,  this  tube  ejected  a  considerable  portion  of  vapor 
bubbles  radially  along  tube  periphery.   This  premature 
departure  of  vapor  bubbles  results  in  a  lower  performance 
by  producing  a  thicker  liquid  film  at  the  hot  surface. 

In  comparison  with  the  results  for  tubes  5,  6  and  7,  the 
curves  representing  tubes  2,  3  and  4  show  that  for  a  given 
fin-tip  gap,  the  performance  has  decreased  with  decrease  in 
fin  density.   This  observation  is  believed  to  be  due  to  the 
presence  of  a  greater  number  of  channels  (for  a  given  length 
of  tube)  in  the  higher-fin-density  tubes.   Again,  for  a  fin 
density  of  0.75  fins/mm,  the  0.25  mm  fin-tip  gap  (tube  3) 
resulted  in  the  best  boiling  performance.   The  same 
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explanations  given  above  should  be  valid  for  the  poorer  per- 
formance of  tubes  2  and  4. 

It  can  be  seen  that  all  tubes  show  smaller  enhancement 
ratios  at  low  wall  superheats,  and  the  above-mentioned 
differences  in  performance  are  more  pronounced.   On  the 
other  hand,  all  the  tubes  show  greater  enhancements  at  larger 
wall  superheats  and  they  seem  to  approach  a  common  enhance- 
ment ratio.   The  increase  in  enhancement  is  believed  to  be 
due  to  the  increased  chimney  effect,  which  helps  to  remove 
vapor  bubbles  from  channels  resulting  in  a  thinner  liquid 
film  at  the  hot  boiling  surface.   When  heat  flux  is  further 
increased,  vapor  can  provide  a  blanket  around  the  hot 
surface,  which  may  eventually  result  in  a  critical  heat 
flux.   It  appears  that  the  critical  heat  flux  is  not  strongly 
dependent  on  the  fin  density  or  fin-tip  gap. 

The  results  shown  in  Figure  28  were  also  replotted  as 
shown  in  Figure  29  to  further  study  the  performance  of  these 
six  tubes.   This  figure  shows  the  variation  in  enhancement 
ratio  obtained  beyond  the  area  increase  (i.e.,  the  ratio  of 
heat-transfer  enhancement/area  ratio)  as  a  function  of  wall 
superheat.   It  is  evident  that  obviously  all  tubes  resulted 
in  enhancement  ratios  beyond  the  area  increase. 

In  summary,  a  fin  density  of  1.02  fins/mm  resulted  in 
performance  superior  to  a  fin  density  of  0.75  fins/mm.   Also, 
in  both  cases,  a  fin-tip  gap  of  0.25  mm  resulted  in  the  best 
performance . 
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V.   CONCLUSIONS 

(1)  Experimental  observations  revealed  a  non-uniform 
temperature  distribution  along  test  section.   The  bubble 
generation  was  seen  to  start  at  the  two  ends  and  migrate  into 
the  center  as  the  heat  flux  was  increased.   No  clear  explana- 
tion is  possible  for  this  temperature  distribution. 

(2)  Results  of  the  present  investigation  showed  very 
good  repeatability  with  the  results  of  Hernandez  obtained 
under  similar  conditions. 

(3)  The  use  of  shrouds  improved  the  boiling  heat- transfer 
coefficient  of  Gewa-T  surface  at  all  practical  heat  fluxes. 

However,  a  considerable  reduction  in  heat-transfer  coefficient 

2 
was  seen  at  very  high  heat  fluxes  (>  100,000  W/m  ) .   Among 

the  four  shrouds  tested,  the  shroud  with  60°  *8.5°  openings 

resulted  in  the  best  performance  verifying  Hernandez's 

results.   This  shroud  resulted  in  a  253  percent  increase  in 

boiling  heat-transfer  coefficient,  over  smooth  tube,  at  a 

wall  superheat  of  1.8  K. 

(4)  The  use  of  0.1  mm-diameter  wires  (2  to  5)  improved 
boiling  heat-transfer  coefficient  at  all  practical  heat 
fluxes.   The  use  of  3  wires  resulted  in  the  best  performance 
with  a  341  percent  increase  in  the  heat-transfer  coefficient 

(over  the  smooth-tube  value)  at  a  wall  superheat  of  1.4  K. 

(5)  Further  improvement  in  boiling  performance  was 
obtained  by  using  the  60°  *8.5°  shroud  with  three  wire  wraps. 


7 


This  combination  resulted  in  a  481  percent  increase  in  boil- 
ing heat-transfer  coefficient  at  a  wall  superheat  of  1.1  K. 

(6)  The  use  of  wires  appears  to  be  superior  to  the  use 
of  more  expensive  shrouds. 

(7)  Among  the  fin-tip  gaps  of  0.15,  0.25  and  0.35  mm 
tested,  the  0.25  mm  gap  resulted  in  the  best  boiling 
performance . 

(8)  The  Gewa-T  tubes  with  a  fin  density  of  1.02  fins/mm 
performed  superior  to  a  fin  density  of  0.75  fins/mm.   At  a 
wall  superheat  of  3.5  K,  the  lower-fin-density  tube  increased 
the  boiling  heat-transfer  coefficient  by  80  percent,  while 
the  higher-fin-density  tube  increased  by  103  percent. 
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VI.   RECOMMENDATIONS 

(1)  Investigate  the  reason  for  the  existence  of  a 
longitudinal  temperature  profile  and  take  appropriate 
steps  to  alleviate  this  problem. 

(2)  Study  the  possibility  of  adding  roughness  in  the 
inter-fin  cavity  of  Gewa-T  tubes.   This  operation  should 
be  performed  before  the  rolling  process  to  obtain  the  T- 
shaped  fins. 

(3)  Repeat  runs  made  during  this  investigation  with 
water  as  the  boiling  medium. 
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APPENDIX  A 
DATA  REDUCTION  PROGRAM 


'000! 

1010! 

i  020! 

1030! 

1040 

1050 

I  060 

1070 

1080 

1090 

1  100 

1110 

I  120 

1130 

1  140 

1  1  50 

1  160 

1170 

I  180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1250 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 


~ILE  "IAME:  DRP2 

DATE:      May  18.  1984 

REVISED:    July  20.  !384 


"4X. 

"6X 
"SX, 


"Select  option:' 
'0  Tak inq  data  or 
'1  Plotting  data'" 

Main 

Plot 


re-orocess inq  previous  data' 


TUBE  QD  (0=0.75  in. 


0  in)". ltd 


BEEP 

PRINTER  IS 

PRINT  USING 

PRINT  USING 

PRINT  USING 

INPUT  Idp 

IF  ldp=0  THEN  CALL 

IF  Idp-1  THEN  CALL 

END 

SUB  Main 

COM  /Cc/  C(7) 

DIM  Emf (9> ,T(9) 

DATA  0. 10086091 .25727.34369.-767345.3295.78025595.31 

DATA  -9247486589.  6. 97683E+1 1  . -2 .661 92E  +  1 3 . 3 . 94073E  +  1 4 

READ  C<*) 

BEEP 

INPUT  "SELECT 

IF  Itd-0  THEN 

D1 - .0159     ! 

D2= . 01 93      !  Uiame 

Di= .0133     ! 

Do=.0190     ! 

L=.0490      ! 

Lu=.0319      !  Length  of  unenhanced  surface  at  the  ends 

ELSE 

D1-.  01786 

D2= .0231 

Di=. 02032 

Do-. 02222 

L=.0503 

Lu=.03175 

END  IF 

Lc=Lu+<Do-Di)/2 

A=PI*(Do" 2-Di  2)/4 

P=PI*Do 

Kcu=385       !  Thermal  conductivity  of  Copper 

PRINTER  IS  701 

CLEAR  709 

BEEP 

INPUT  "ENTER  MONTH.  DATE  AND  TIME  ( MM : DD : HH :MM : SS>" . DateS 

OUTPUT  709 :"TD": DateS 

OUTPUT  709:"TD" 

ENTER  709: DateS 

Month,  date  and  time 


Diameter  at  thermocouple  positions 
Diameter  of  test  section  to  the  base 
Inside  diameter  of  unenhanced  ends 
Outside  diameter  of  unenhanced  ends 
Length  of  enhanced  surface 
Length  of  unennanced  surface  at  the 


PRINT 

PRINT 

PRINT 

BEEP 

INPUT 

PRINT 

BEEP 

INPUT 


USING  "10X.""NOTE:  Program  name 


' : DateS 
BOIL 


"ENTER  DISK 
USING  "1GX." 


NUMBER", Dn 
"Disk  number 


"ENTER  INPUT  MODE  (0= 


""  .ZZ":Dn 
3054A.  I -FILE)", In 


Im  = 

IF 


Im+1 


THEN 


BEEP 

INPUT  "GIVE 

PRINT  USING 


A  NAME 
"16X." 


P0R 
'New 


THE  RAW  DATA 
file  name : 


FILE" ,D_f LieS 

.  14A" :D  f ileS 
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1600 
!6!0 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
!  700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
!300 
1810 
1320 
1830 
1840 
1850 
1860 
1870 
1380 
1390 
1900 
1910 
1920 
1930 
1940 
1950 
i960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2'60 
21  70 
2!  80 
2190 


BEEP 

INPUT  "ENTER  FLUID  ( 0=FREON- 113 . 1 -WATER)" .Cf 

Cf-Cf+1 

CREATE  BDAT  D_FileS,20 

ASSIGN  $File  TO  D_f ile£ 

OUTPUT  ?iFile:Date5,Cf 

ELSE 

BEEP 

INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" .D_f i leS 

BEEP 

Ihc-0 

INPUT  "ENTER  1  IF  USING  FILE  NUMBER  <  25" .Ihc 

PRINT  USING  "16X.,,M0ld  file  name:  ""  .  1  4A"  :  D_f  l  leS 

BEEP 

INPUT  "ENTER  THE  NUMBER  OF  RUNS  STORED"  .Nrtin 

ASSIGN  «File  TO  D_file$ 

ENTER  aFile:DoidS,Cf 

PRINT  USING  "16X,""This  output  for  file  created  at  ""  .  1 4A" : DoIdS 

END  IF 

BEEP 

INPUT  "GIVE  A  NAME  FOR  PLOT  FILE" .P_f i leS 

CREATE  BDAT  P_fiie$.5 

ASSIGN  aPlot  TO  P_fileS 

IF  Cf-1  THEN  PRINT  USING  "16X .""Fluid 

IF  Cf-2  THEN  PRINT  USING  "16X. ""Fluid 

BEEP 

INPUT  "ENTER  OUTPUT  VERSION  < t -SHORT. 0-LQNG)" . Iov 

IF  lov=0  THEM  Iov=2 

BEEP 

INPUT  "ENTER  NUMBER  OF  DEFFECTIVE  TCS  <Q-DEFAULT>' 

IF  ldtc>0  THEN 

BEEP 

IF  Idtc=1  THEM 

INPUT  "ENTER  DEFLECTIVE  TC  LOCATION" .Ldtcl 

PRINT  USING  "16X.""TC  is  deffective  at  location  "' 

END  IF 

IF  Idtc-2  THEN 

INPUT  "ENTER  DEFFECTIVE  TC  LOCATIONS" , Ldtc 1 ,Ldtc2 

PRINT  USING  "16X.""TC  are  deffective  at  locations 

END  IF 

END  IF 

IF  Im-f  THEN 

BEEP 

INPUT  "ENTER  BAROMETRIC  PRESSURE  (inHg)".Patm 

Patm=Patm+2*1 500/1 3500 

OUTPUT    2>Fiie:Patm 

ELSE 

ENTER    »)File:Patm 

IF  Ihc-1  THEN  Patm=Patm+2*1500/13500 

END  IF 

Patm=Patm*25 . 4* 1 4 . 696/760 

Tsat=FNTsat<Cf .Patrn) 

PRINT  USING  "16X. ""Saturation  temperature  =  "" . 3D . DD . ""  < Deg  C)' 

J-1 

Sx=0 

Sy  =  0 

Sxs  =  0 

Sxy  =  0 

BEEP 

INPUT  "HIT  CONTINUE  WHEN  READY".  0k 


Freon-1 1 3' 
Water 


Idtc 


D":Ldtc1 


D.4X.D":Ldtc1 .Ldtc2 


;  Tsat 
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2200  Repeat:  ! 

2210  IF  Im-I  'HEM 

2220  OUTPUT  709:"AR  AF29  AL33  VR5" 

2230  FOR  r=0  TQ  9 

2240  OUTPUT  709: "AS  SA" 

2250  Suin=0 

2260  FOR  J i=1  TO  20 

2270  ENTER  709:E 

2280  Smn  =  Som+E 

2290  NEXT  Ji 

2300  Emf<I>=Sum/20 

2310  NEXT  I 

2320  BEEP 

2330  INPUT  "ENTER  HEATER  VOLTAGE.  Vh  (V),,.Vh 

2340  BEEP 

2350  INPUT  "ENTER  RESISTOR  VOLTAGE.  Ur  (V>".Vr 

2360  ELSE 

2370  ENTER  aFile:Emf <*> .Vr.Vh 

2380  END  IF 

2390  Tua-0 

2400  FOR  1=0  TO  9 

2410  IF  ldtc>0  THEN 

2420  IF  (I+1)-Ldtc1  OR  <I+1)-Ldtc2  THEM 

2430  T(I)=99.99 

2440  GOTO  2490 

2450  END  IF 

2460  END  IF 

2470  T(I)=FNTvsv(Emt  <D  ) 

2480  IF  I<8  THEN  Tua=Twa+T(  I ) 

2490  NEXT  I 

2500  Twa=Tuia/<8-Idtc) 

2510  Tld=T(9) 

2520  Tv=T<8> 

2530  Rho=FMRho<Cf .Tsat) 

2540  Mu=FNMu(Cf .Tsat) 

2550  K=FNK(Cf .Tsat) 

2560  Cp=FNCp(Cf .Tsat) 

2570  8eta=FNBeta(Cf .Tsat) 

2580  Ni=Mu7Rho 

2590  Alpha=K/(Rho*Cp) 

2600  Pr=Ni/Alpha 

2610  Amp=Vr/2.031 

2620  Q=Vh*AmD 

2630  Tu-Tua-Q-LQG<D2/D1 > /<2*PI*Kcu~L ) 

2640  Thetab=Tw-Tsat 

2650  IF  ThetaD<0  THEN 

2660  BEEP 

2670  OISP  "TWALL<TSAT  (HIT  CONTINUE  WHEN  READY)" 

2680  PAUSE 

2690  GOTO  2210 

2700  END  IF 

2710  Hbar=190 

2720  Fe=<Hbar*P/(Kcu.*A)  >  .5*Lc 

2730  Tanh=FNTanh(Fe> 

2740  Theta=Thetab-Tann/Fe 

2750  Xx  =  <9.31*Beta*Thetab*Do   3*Tanh/<Fe*Ni*fHpha>  )     .166667 

2760  Yy=(1+(  .559/Pr)    (9/16) )' (8/27) 

2770  Hbarc=K/Do*<  .6+.387»Xx/Yy)  2 

2780  IF  ABS< < Hbar-Hbarc) /Hbarc ) > . 00 1  THEM 

2790  Hbar=(Hbar+Hbarc;*.5 

2800  GOTO  2720 


2810  END  IF 

2320  Ql  =  (Hbar*P*x.cu*A)  . 5*Thetab*Tann 

2830  Qc=Q-2*Ql 

2840  fls=PI*D2*L 

2850  Qdp=Qc/As 

2360  IF  Iov-I  THEN  PRINT  USING  '"!  OX  .  DD  .  2X  .  3D  .  DD  .  2X  .2<  Z  .  4DE  .  2X>  .  DD  .  DD"  :  J  .  Twa.Hba 
r .Qdp. Thetab 

.2370  IF  Iov=2  THEN 

2380  PRINT  USING  "1  OX. ""Data  set  number  =  ""  .  DD"  :  J 

2890  PRINT  USING  "1  OX  .  ,,"Posi  t  ion  numoer  :    1      2     3     4     5     6 

1  8 

2900  PRINT  USING  "1  OX . ""Temperature     :  ""  .  3( DD . DD . I  X) " : T < 0 ) . T< 1) . T( 2 ) , T( 3 ) . T 
(4)  .  T<5> ,T<6> .T(7> 

2310  PRINT  USING  "1  OX  .  ""Ti igd    Tvapr    Twa      Hbar        Qdp     Thetab 

2920  PRINT  USING  "3X ,3(3D . DD .2X) . 2( Z . 4DE .2X > . 2D . DD" : T Id . Tv . Twa. Hbar , Qdp . ThetaD 

2930  END  IF 

2940  IF  Im-1  THEM 

2950  BEEP 

2960  INPUT  "OK  TO  STORE  THIS  DATA  SET  (  1  -Y  .  0=N>?" .Ok 

2970  END  IF 

2980  IF  Ok-1  OR  Im-2  THEN  J  =  J  +  1 

2990  IF  Ok-1  AND  Im-1  THEN  OUTPUT  @Fi  le  :Emf  <  *  >  .  Vr  .  Vh 

3000  IF  Im=2  OR  Ok=1  THEN  OUTPUT  «P lot :Odp . Thetab 

3010  IF  Im-1  THEN 

3020  BEEP 

3030  INPUT  "WILL  THERE  BE  ANOTHER  RUN  ( I -Y . 0-N) 7" ,Go_jn 

3040  Nrun=J 

3U50  IF  Go.onOI  THEN  3110 

3060  WAIT  300 

30  70  IF  Go_on=!  THEN  Repeat 

3080  ELSE 

3090  IF  J<Nrun+1  THEN  Repeat 

3100  END  IF 

3110  IF  Im-1  THEM 

3120  BEEP 

3130  PRINT 

3140  PRINT  USING  "t  OX .  ""NOTE:  "••.ZZ,'■•,  data  runs  were  stored  in  file  ,,,,,!0A";J- 

1 .D_f ileS 

3'50  END  IF 

3160  BEEP 

3  t  70  PRINT 

3180  PRINT  USING  "1  OX . ""NOTE :  "".ZZ.""  X-Y  pairs  were  stored  in  plot  data  file 

"", 10A":J-1 ,P_f ileS 

3190  ASSIGN  iFile  TO  * 

3200  ASSIGN  wPlot  TO  * 

3210  BEEP 

3220  INPUT  "LIKE  TO  PLOT  DATA  ( 1 -Y , 0»N>?" .Ok 

3230  IF  Ok-1  THEM 

3240  CALL  Plot 

3250  END  IF 

3260  SUBEND 

3270  DEF  FNMu<Cf . T> 

3230  IF  Cf-I  THEM 

3290!  RANGE:  28  TO  50  DEG  C 

3300  ttu=8.3629813E-4-T*< 1 . ! 094609E-5-T«5 . 566829E-8) 

3310  END  IF 

3320  IF  Cf-2  THEN 

3330  A=247.3/<T+133.15> 

3340  Mu=2.4E-5*l0  A 

3350  END  IF 
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3360  RETURN  liu 

3370  FNEND 

3380  DEF  FNCp(Cf.T) 

3390  IF  Cf-1  THEN 

3400!  RANGE:  0  TO  30  DEG  C 

3410  Cp=9.2507273E-!+T*<9.3400433E-4+1 . 7207792E-6*T) 

3420  END  IF 

3430  IF  Cf-2  THEN 

3440  Cp=4.21120858-T*<2.26826E-3-T*(4.42361E-5+2.71428E-7*T) ) 

3450  END  IF 

3460  RETURN  Cp*1000 

3470  FNEND 

3480  DEF  FNRho(Cf.T) 

3490  IF  Cf-1  THEN 

3500!  RANGE:  30  TO  80  DEG  C 

3510  Ro=1  .6207479E+3-TM2.2186346  +  T*2.3578291E-3> 

3520  END  IF 

3530  IF  Cf-2  THEN 

3540  Ro=999.52946  +  T*(  . 0 !269-t*<5. 48251 3E-3-T-1  .234M7E-5)) 

3550  END  IF 

3560  RETURN  Ro 

3570  FNEND 

3580  DEF  FNPr(Cf.T) 

3590  Pr=FNCp<Cf ,T)*FNMu(Cf . T ) /FNK ( Cf . T) 

3600  RETURN  Pr 

3610  FNEND 

3620  DEF  FNK(Cf.T) 

3630  IF  Cf-1  THEN 

3640!  RANGE:  30  TO  90  DEG  C 

3650  K  =  8.2095238E-2-TM2.2214286E-4  +  T»2.3809524E-8> 

3660  END  IF 

3670  IF  Cf-2  THEN 

3680  X=(T+273. 15>/273.15 

3690  K=-.92247+X*<2.8395-X*< I  ,3007-X*< . 52577- . 07344»X  ))  ) 

3700  •  END  IF 

3710  RETURN  K 

3720  FNEND 

3730  DEF  FNTanh(X) 

3740  P=EXP(X) 

3750  Q=EXP(-X> 

3760  Tanh=<P-0)/(P+Q) 

3770  RETURN  Tanh 

3780  FNEND 

3790  DEF  FNTvsv(V) 

3800  COM  /Cc/  C(7) 

3310  T=C<0) 

3820  FOR  1=1  TO  7 

2830  T-T+C(I>»V"I 

3840  NEXT  I 

3850  T=T+8.34643E-l-T*<2.4547896E-2-T»! .8722583E-4) 

3860  RETURN  T 

3370  FNEND 

3880  DEF  FNBeta(Cf.T) 

3890  Rop=FNRho(Cf ,T+. 1 > 

3900  Rom=FNRho<Cf  .T-. 1 ) 

3310  Beta=-2/<Rop+Rom)*<Rop-Rom>/  .2 

3920  RETURN  Beta 

3930  FNEND 

3940  DEF    FNTsaUCf .P> 

3950  Tl-10 

3960  Th-110 


3970 
3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4061 

4070 

4080 

4081 

4082 

4083 

4085 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4I60 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 


2635*LGT(T>+2.0539E-3»T) 


Ta=<Tl+Th>*.5 

Pc=FNPsat<Cf  .  Ta) 

IF  ABS< (P-Pc)/P)>.0001  THEf 

IF  Pc>P  THEN  Th=Ta 

IF  Pc<P  rHEN  Tl=Ta 

GOTO  3970 

END  IF 

RETURN  Ta 

FNEND 

DEF  FNPsat(Cf . Tc) 

IF  Cf-1  THEN 

T=Tc*! .8+32+459.6 

P=I0"<33. 0655-4330. 98/T-9 

END  IF 

IF  Cr=2  THEN 

P=Pvstu(Tc) 

END  IF 

RETURN  P 

FNEND 

SUB  Plot 

DIM  C(9) 

INTEGER  Ii 

PRINTER  IS  1 

Idv-1 

BEEP 

INPUT  "LIKE  DEFAULT  VALUES  FOR  PLOT  <!=Y.0=N> 

BEEP 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

INPUT 

IF  0Po=3 

BEEP 

INPUT  "SELECT 

END  IF 

PRINTER  IS  705 

IF  IdvOI  THEN 

BEEP 

'ENTER 


Idv 


USING 
USING 
USING 
USING 
USING 
Opo 

THEN 


"4X, 

"4X 

"4X, 

"4X 

"4X, 


"Select  Option: 

"0   q  versus  delta-T' 

h   versus   delta-T 

h    versus    q 

h-ratio    versus    delta-T' 


'1 


TUBE  DIAMETER  <0=.75 


0  IN> ".ltd 


INPUT 

BEEP 

INPUT  "ENTER 

BEEP 

INPUT  "ENTER 

BEEP 

INPUT  "ENTER 

ELSE 

Cy  =  3 

IF  0do=0  THEM 

Cx-2 

Xmi n= 1 

Ymin=!000 

END  IF' 

IF  Odo-I  THEN 

Cx=2 

Xm l n= ' 

Ynin=1 00 

END  IF 

IF  0po=2  THEM 

Cx=3 


NUMBER 
NUMBER 


OF 

OF 


CYCLES 
CYCLES 


FOR 
FOR 


X-AXIS". Cx 


Y-AXIS". Cy 
MIN    X-VALUE    (MULTIPLE   OF    IO)",Xmin 
MIN    Y-VALUE    (MULTIPLE   OF    IO>",Ymin 


86 


4530 

Xmin=l00G 

-4540 

Ymin='00 

4550 

END  IF 

4b60 

IF  dpo-3  THEN 

4570 

Cx  =  2 

4580 

Cy  =  3 

4590 

Xirun=1 

4600 

Ymin= 1 

4610 

END  IF 

4620 

END  IF 

4630 

BEEP 

4640 

PRINT  "IN:SP1 :IP  2300. 18 

00,3300 

.6800 

4650 

PRINT  "SC  0  . 1 00  .  0  . 1 00 : TL 

2.0;" 

4660 

Sfx=100/Cx 

4670 

Sfy=100/Cy 

4680 

PRINT  "PU  0.0  PD" 

4690 

Nn  =  9 

4700 

FOR  1=1  TO  Cx+1 

4710 

Xat=Xmin*10' ( 1-1 > 

4720 

IF  r-Cx+1  THEN  Nn-1 

4730 

FOR  J=1  TO  Nn 

4740 

IF  J-1  THEN  PRINT  "TL  2 

0" 

4750 

IF  J=2  THEN  PRINT  "TL  1 

0" 

4760 

Xa=Xat*J 

4770 

X=LGT(Xa/Xmin>*Sf x 

4780 

PRINT  "PA":X.".0:  XT:" 

4790 

NEXT  J 

4800 

NEXT  I 

43  10 

PRINT  "PA  I00.0:PU:" 

4820 

PRINT  "PU  PA  0.0  PD" 

4830 

Nn  =  9 

4840 

FOR  1=1  TO  Cy+1 

4850 

Yat  =  Ymin*10'  <I-1 ) 

4860 

IF  I=Cy+!  THEN  Nn=1 

4370 

FOR  J-1  TO  Nn 

*380 

IF  J-1  THEN  PRINT  "TL  2 

0" 

4890 

IF  J-2  THEN  PRINT  "TL  1 

0" 

4900 

Ya=Yat*J 

4910 

Y  =  LGT(Ya/Ynun>»Sfv 

4920 

PRINT  "PA  0.":Y."YT" 

4930 

NEXT  J 

4940 

NEXT  I 

4950 

PRINT  "PA  0,  100  TL  0  2" 

4960 

Nn  =  9 

4970 

FOR  1=1  TO  Cx+1 

4980 

Xat  =  Xmin*10'  ( 1-1 ) 

4990 

IF  I-Cx+t  THEN  Nn=1 

5000 

FOR  J-1  TO  Nn 

5010 

IF  J-1  THEM  PRINT  "TL  0 

->  *• 

5020 

IF  J>1  THEN  PRINT  "TL  0 

1" 

5030 

Xa=Xat*J 

5040 

X  =  LGT(Xa/Xnun>*Sfx 

5050 

PRINT  "PA";X.", 100:  XT" 

5060 

NEXT  J 

5070 

NEXT  I 

5080 

PRINT  "PA  100.  100  PU  PA 

100.0  P 

D" 

5090 

Nn  =  9 

5'00 

FOR  1=1  TO  Cy+! 

5110 

Yat  =  Ynun~10  ( 1-1 > 

5'20 

It-  I=Cy+!  THEN  Nn-! 

5130 

FOR  J-1  TO  Nn 
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5140  IF  J- I  THEN  PRINT  "TL  Q  2" 

5150  IF  J>1  THEN  PRINT  "TL  0  1" 

5160  Ya-Yat*J 

5170  Y=LGT(Ya/Ymin)*Sfy 

5130  PRINT  "PD  PA  I  00 ." . Y ."YT" 

5190  NEXT  J 

5200  NEXT  I 

5210  PRINT  "PA  100.  !U0  PIT 

5220  PRINT  "PA  0.-2  SR  1  .5,2" 

5230  Ii=LGT(Xmin) 

5240  FOR  1=1  TO  Cx  +  1 

5250  Xa=Xmin*10  (1-1 > 

5260  X=LGT(Xa/Xmin>*Sfx 

5270  PRINT  ,,PA":X.".0:" 

5290  IF  li>  =  0  THEN  PRINT  "CP  -2 , -2 :LB1 0 :PR  -2  .2  ;LB"  ;  1 1  ;  "" 

5290  IF  li<0  THEN  PRINT  "CP  -2  .  -2 :LB1 0 :PR  0 ,2;LB" : Ii :"" 

5300  Ii-Ii+1 

5310  NEXT  I 

5320  PRINT  "PU  PA  0.0" 

5330  Ii-LGT(Ymin) 

5340  Y10-I0 

5350  FOR  1=1  TO  Cy  +  1 

5360  Ya=Ymin*10  (1-1) 

5370  Y=LGT(Ya/Ymin)*Sfy 

5380  PRINT  "PA  0.":Y,"" 

5390  PRINT  "CP  -4  .  -  .25  :  LB1 0  :PR  -2 .2 ; LB"  :  1 1 : "" 

5400  Ii-Ii+1 

5410  NEXT  I 

5420  IF  IdvOI  THEN 

5430  BEEP 

5440  INPUT  "ENTER  X-LABEL" . X laDeiS 

5450  BEEP 

5460  INPUT  "ENTER  Y-LABEL" . Y laoelS 

5470  ELSE 

5480  IF   Opo-0    THEN 

5490  XlabelS»"Tuj-Tsat    (K>" 

5500  YlabelS="q    <W/m   2)" 

5510  END    IF 

5520  IF  Dpo=1  THEN 

5530  XlabelS="Tu-Tsat  < K ) " 

5540  YlabelS="h  <W/m  2.K)" 

5550  END  IF 

5560  IF  0po=2  THEN 

5570  XlabelS="a  (W/m  2>" 

5580  YlabelS="h  (W/m  2.K)" 

5590  END  IF 

5600  IF  Gpo=3  THEN 

5610  XlabelS="Tu-Tsat  <K>" 

5620  Y Labe lS="h(enh)/h( smooth)" 

5630  END  IF 

5640  END  IF 

5650  PRINT  "SR  l.5.2;PU  pA  50.-16  CP" : -LEN(XlabelS>/2: "0 :LB" : XlabelS : ' 

5660  PRINT  "PA  -14.50  CP  0 . " : -LEN< Y labe IS) /2*5/6 : "DI  0 . 1  :LB" : Y labe IS : ' 

5670  PRINT  "CP  0.0  DI" 

5680  Repeat: ! 

5690  BEEP 

5700  INPUT  "WANT  TQ  PLOT  DATA  FROM  A  FILE  (  I -Y  ,0-N)?" .Ok 

5710   IF  Q|c-1  THEM 

5720  BEEP 

5730   INPUT  "ENTER  THE  NAME  OF  "HE  DATA  FILE" .D_f i leS 

5740  ASSIGN  ©File  TO  D  f ileS 


5750 
5760 
5770 
5780 
5790 
5800 
5810 
5320 
5830 
5840 
5850 
5860 
5870 
5380 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
6060 
5070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
6150 
6160 
6170 
6130 
6190 
6200 
6210 
6220 
6230 
6240 
5250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 


BEEP 
BEEP 
INPUT 
BEEP 
INPUT 
BEEP 

PRINTER  IS 
PRINT  USING 


•ENTER 
•ENTER 


THE  BEGINNING  RUN  NUMBER" . Md 

THE  NUMBER  OF  X-Y  PAIRS  STORED" .Npairs 


USING 
USING 
USING 
USING 
USING 
Sym 


"4X 
"4X 
"4X 
"4X 
"4X 
"4X 


elect  a  symbol  : '" 

Star      2   PI  as  sign"' 
Circle    4   Square""" 

Rombus 

Right-side-up  triangle' 
Up-side-doun  triangle'" 


PRINT 

PRINT 

PRINT 

PRINT  USING  "4X.""6 

PRINT  USING  "4X,""7 

INPUT 

PRINTER  IS  705 

PRINT  "PU  DI" 

IF  Sym=1  THEN  PRINT  "SM»" 

IF  Sym=2  THEN  PRINT  "SM+" 

IF  Sym=3  THEN  PRINT  "SMo" 

IF  Md>1  THEN 

FOR  1=1  TO  <Md-1 ) 

ENTER  «File;Ya.Xa 

NEXT  I 

END  IF 

FOR  1=1  TO  Npairs 

ENTER  §»File:Ya.Xa 

IF  0po=1  THEN  Ya=Ya/Xa 

IF  0po=2  THEN 

Q  =  Ya 

Ya=Ya/Xa 

Xa=0 

END  IF 

IF  0po=3  THEN  Ya=Ya/FNHsmooth(Xa 

X=LGT(Xa/Xmin)*Sfx 

Y=LGT(Ya/Ymin)*Sfy 

IF  Sym>3  THEN  PRINT  "SM" 

IF  Sym<4  THEN  PRINT  "SR  1.4,2.4" 

PRINT  "PA".X.Y."" 

IF  Sym>3  THEN  PRINT  "SR  1 

PRINT 

PRINT 

PRINT 

PRINT 


ltd) 


IF 
IF 
IF 
IF 


THEM 
THEN 
THEN 
THEN 


"SR 

"UC2 
"UC3 
"IJCO 
"UC0 


4.99 
0.99 


6" 

0.-8.-4.0.0.3.4.0 
.-3,-6,-3.6.3.6.3.' 
99.3.-8.-6.0,3.8:" 


5.3.99.-3.3.6.0.-3.-8:" 


Sym  =  4 
Sym  =  5 
Sym  =  6 
Sym=7 

NEXT    I 

BEEP 

ASSIGN  &File  TO  * 

GOTO  5690 

END  IF 

PRINT  "PU  SM" 

BEEP 

INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  H=Y.0=N) 

IF  Go_on=1  THEN 

BEEP 

INPUT  "ENTER  LOWER  AND  UPPER 

FOR  Xx  =  0  TO  Cx  STEP  Cx/200 

Xa=Xmin*l0  Xx 

IF  Xa<Xll  OR  Xa>Xlu  THEN  6330 

Ya=FNPoly(Xa) 

Y=LGT< Ya/Ymin)»Sfy 

X=LGT(Xa/Xmin)*Sfx 

IF    Y<0    THEN    Y=0 


,  bo    on 


X-LIMITS".Xll,Xlu 
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536(J  ir  Y >  i  0 U  "HEN  GOTO  6330 

6370  PRINT  "PA".X.Y."PD" 

6330  NEXT  Xx 

6390  END  IF 

6400  PRINT  "PU  PA  0.0  3P0" 

6410  SUBEND 

6420  DEF  FNHsmoorhCX . ltd) 

6430  Hs=FNPoly<X>/X 

6440  IF  Itd=1  THEN  Hs=Hs*. 83347 

6450  RETURN  Hs 

6460  EMEND 

6470  DEF  FNPoly(X) 

6480  Poly=-4.41237I8E+2-X»<6.8123917E+2-X»3.7416863E+2> 

6490  RETURN  Poly 

6500  FNEND 

6510  DEF  FNPvst(Tsteajn) 

6515  DIM  K<8> 

6520  DATA  -7 . 691 234564 . -26 . 08023696 .-! 68 . 1706546.64.23285504.-1 18.9646225 

6525  DATA  4 . 1 671 1 732  .20 .9750676 . 1  ,E9 ,6 

6530  READ  K<») 

6535  T=(Tsteam+273.15)/647.3 

6540  Sum=0 

6545  FOR  N=0  TO  4 

6550  Sum  =  Sum+K<N>*<  l-T)  (N+!  ) 

6555  NEXT  M 

6560  Br=Sum/<T*<  1+K<5)*<1-T)-H<<6>*<1-T>  "2>  >-<  1  -T>/(K<7>*<  ?-T)  "2+K(8>) 

6565  Pr=EXP<Br> 

6570  P=22120000*Pr 

6575  RETURN   P 

6580  FNEND 
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APPENDIX  B 
SAMPLE  CALCULATIONS 

Data  set  number  fifteen  of  file  GT6H29  corresponding  to 
tube  number  6  was  chosen  for  the  sample  calculation. 

A.   TEST-SECTION  DIMENSIONS 


O.D. 

= 

25.30 

(mm) 

HF 

= 

1.05 

(mm) 

fins/mm 

= 

1.02 

inter- 

-fin 

gap 

= 

0.25 

(mm) 

D2 

= 

23.10 

(mm) 

Dl 

= 

17.86 

(mm) 

D 
o 

= 

22.22 

(mm) 

D. 

l 

= 

20.32 

(mm) 

Lu 

= 

31.75 

(mm) 

L 

= 

50.80 

(mm) 

B.   MEASURED  PARAMETERS 

V  =   70.0  (volts) 


VRES  =    4.8  (volts) 

R  =2.0  31  (ohms) 

kc  =  383  (W/m.K) 

T1  =      54.23  (°C) 

T2  =   54.29  (°C) 

T3  =   53.91  (°C) 
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T4  =  53.65  (°C) 

T5  =  53.40  (°C) 

Tc  =  53.31  (°C) 

T?  =  53.22  (°C) 

Tg  =  53.54  (°C) 

TSAT  '   =  47.51  CC) 


C.  PROPERTIES  OF  FREON-113  AT  SATURATION  TEMPERATURE 
U  =    4.948  x 10"4  (kg/m.s) 

p  =   1510.0  (kg/m3) 

Cp  =   0.973  (kJ/kg.K) 

k  =   0.0715  (W/m.K) 

D.  HEAT-FLUX  CALCULATION 

The  heat  transfer  rate  from  the  cartridge  heater: 


QH   =   VI  (W)     and     I   =   VR£S/R  (amps) 


I   =   4.8/2.031   =   2.363  (amps) 


Qu      =      70.0x2.363      =      165.410     (WO 
H 


The  cross-sectional  area  of  machined  ends 


An   =   n/4.0(D2  -  D2)  (m2) 
C  o     i 

A        =      tt/4.0  (0.022222    -    0.020322)       =      63.481  xio"6     (m2 
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The    tube    outside   wall    perimeter: 


P      =      uD         =      3.1416x0.02222      =       6.981x10    2     (m) 
o 


The  corrected  length  of  ends 


Ln      =      Lr7   +    t/2.0  t      =      D      -    D. 

C  U  o  l 

L        =      0.03175    +    (0.02222    -    0.02032)/2.0      =      32.7  x 10"3     (m) 


Other  thermophysical  properties: 


v      =      y/p      =      4.948x10    4/1510.0      =       3.277x10    7     (m2/s; 


a      =      k/pCp      =  0.0715/(1510.0x973)       =      4.866x10    8     (m2/s; 

1    A_p  1      ,1509.77  -  1510.26, 

y                   p    AT  15101       47.61    -    47.41    ' 

=  0.00162     (1/K) 


Pr      =      v/a      =      3.277x10    7/4. 866x10    8      =       6.734 


Temperature   at    the  base   of    the    T-shaped    fins 


ln(D   /D    ) 
Tn       =       T_„_    -    QT 


B  AVG         *H         2  tt  L    k, 


TAVG       =        (     I       TN)/8-°       =       53'69     (°C: 
n=l 


T  53    69    -    165    41    ln (23- 1/17 . 86)  5_    _.        ,0    . 

B  -LDD.ix    2ttx0. 0508x383  DJ,J4       K    U' 
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B       =       TB"TSAT       =       53'342    "    47'51       =       5'832     (R) 


Average  heat-transfer  coefficient  at  the  unenhanced 
ends  : 


^DoeBTH(k^)1/2Lc| 


i  1/6 


vaL_,  ( 


C    C 

hP    , 1/2 


h      =      ^-<  0.60  +0.387 


c  kcAc 


[l+(0.559/Pr)9/16]8/27 


h      =      187.49     (W/m    .K) 


Heat- transfer  rate  through  unenhanced  ends 


Jk 


QF      =        VhPkcAc     9B    Tanh   mLc 


m      = 


hP    ,1/2  .187.489   x6.981  x 10 


-2       1/2 


kcV 


0.0715  x  63.481  x  10 


-6 


=   1.698  x  10" 


QF   =   3.29  (W) 


QL0SS   =   2QF   =   6-579  (W: 


Heat  flux  through  Gewa-T  surface 


Q   "   QH  "  QLOSS 
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Q   =   165.41  -  6.579   =   158.83  (W) 


q   =   Q/Ag 


A        =      tt  x23.1  x  10    3  x50.8  x  10    3      =      3.687   x  10    3     (m2 


=      43.978  x 103     (W/m2 


3.687  x 10  3 


The  following  are  the  results  obtained  from  the  computer 
by  running  the  data  reduction  program  (DRP)  (see  Appendix  A) 


q   =   44.466   10  (W/m2) 


>B   =   5.76  (K) 


h"  =   187.39  (W/m2.K) 


The  small  difference  in  results  obey  the  inaccuracy  of 
hand  calculations  with  the  use  of  about  three  significant 
digits  only. 
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APPENDIX  C 
UNCERTAINTY  ANALYSIS 

The  same  data  set  (set  No.  15  of  file  GT6H29)  that  was 
used  for  the  sample  calculation  was  chosen  for  the  uncertainty 
analysis;  therefore  dimensions  of  test  section,  measured  and 
calculated  parameters  found  in  sample  calculation  were  used 
in  this  analysis.   All  uncertainties  are  presented  as  a 
percentage  of  the  calculated  parameter. 

A.   UNCERTAINTY  IN  SOURCE  HEAT- TRANSFER  RATE 


Q   =   VI  (W) ,       I   =   VRES/R  (amp) 


6V   =   SV^g   =   0.05  (volts) 


where 


6   =   uncertainty  in  measurement  and  calculation 


61      =      5VT5r,c/R      =      0.05/2.031      =      0.025     (amp) 


6QH  ...    2  2    1/2 

"or     =      K-)         (— )    ] 
H  VI 


|M     .      ((«-,^    («|,V/2      =      1. OS   percent 
H 
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B.   UNCERTAINTY  IN  SURFACE  AREA 


A^      =      ttD^L  6D0      =      6T  0.1    (mm) 

D  Z  Z  Li 


5AB  0     1       2  0     1      2     1/2 

aJ    ■      "&T'      +    <M»    ]  =     °-48  Percent 


C.   UNCERTAINTY  IN  AT 


AT   =   T   -  T  6T    =   5T      =   0.15  (°C) 

ai       W     SAT         W        SAT  V   ; 


AT  AT 


5AT      r,0.15N2    ,   0.15,2 


AT"   =  [i67T8)       +    (_  6TT8)  ]      =3.43  percent 


D.   UNCERTAINTY  IN  HEAT  FLUX 


QH  *  2QF 

q    =    ■ 


*B 


Ra                      6QH    2      2SQF      2  SAB    2  1/2 

23,      =   [( 5 )z    +     i EL)''  +  (_  — 2.)zl 


In  this  case,  Qu  =  165.41  (W)   and   Q_  =  3.29  (W) 
n  r 


97 


Therefore,  Q   =  2.0  percent  of  Q  .   Assuming  the  same 

proportion  in  the  uncertainty  for  Q  : 

F 


6Q„      =      0.02    6QU      =      0.02x1.753      =      0.035    (W! 
r  n 


5q      _               1.753    2                 2  x  0.035,2  nn48>2l1/2 

T      "       L(158.83}       +    ( 158.83    )       +    (-°-0048)     J 


—     =      1.26    percent 


E.   UNCERTAINTY  IN  BOILING  HEAT-TRANSFER  COEFFICIENT 


h      = 

q 

AT 

h 

-       [(6a,2    +     (.    5AT      1/2 

q                 v         AT  ;  J 

—  1/2 

—  =   [(0.0126)2  +  (-0.0343)2]      =   3.65  percent 
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